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Abstract 

Some of the properties of the low temperature vortex-glass phase of randomly- 
pinned flux lines in 1+1 dimensions are studied. The flux arrays are found to 
be sensitive to small changes in external parameters such as the magnetic field 
or temperature. These effects are captured by the variations in the magnetic 
response and noise, which have universal statistics and should provide an 
unambiguous signature of the glass phase. 
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Flux lines in a clean Type-II superconductor form an Abrikosov lattice at low temper- 
atures 0. However, the flux lattice is destroyed by random microscopic impurities in the 
material |2j . Recently, it has been suggested that the disordered flux array may form a new 
thermodynamic phase at low temperatures, called the "vortex glass" phase, in which flux 
lines are collectively pinned by the impurities |3|-|5| . Although glass-like behavior has been 
reported experimentally ||, a quantitative theoretical description of the vortex-glass phase 
is still lacking except in the special case of flux lines confined to a plane (1+1 dimensions). 
As first shown in Ref. ||, such a 1 + 1 dimensional flux array undergoes a phase transition 
at a finite temperature T g . Below T g , the flux array is pinned by the random impurities 
and forms a glass phase. However, the properties of the glass phase have not yet been elu- 
cidated, and a number of contradictory results exist in the literature |7|-|13|. In this paper, 



we analyze the vortex glass phase using the renormalization-group method of Cardy and 
Ostlund |14| |. We find the glass phase to be characterized by anomalous variations in the 



magnetic responses of the flux array, and extreme sensitivity to small changes in the ap- 
plied field, impurity potential, and temperature. Such glassy behavior has been previously 
conjectured for spin glasses [jOJ and one flux line [p)|-|T%[], and is expected to be generic to 



a wide class of randomness dominated phases. However, the 1+1 dimensional flux array is 
one of the very few systems where analytic results can be obtained. 

We consider an array of flux lines confined to the (x, z)-plane, with an applied field 
H = H z z and repulsive interactions which we model by linear elasticity |1^,^]. Impurites 



yield a random potential V(x,z). Labeling the transverse displacement of the n-th line by 
r n (z) — (n — (fi n (z) /2n) / p, where p ~ H z is the average line density, we can describe the 
large scale fluctuations of the flux array by the Hamiltonian [j3|,[Fysl,[20f 

H = J dxdz || [(d x (f)) 2 + (d z (f)) 2 ] - V(x, z)d x <p - AnpVcos [2irpx + 0] J , (1) 

where <ft(x, z) is the coarse-grained displacement field. In Eq. ([[]), the x— and z— di- 
mensions have been rescaled to make the quadratic part isotropic. The elastic coefficient 
k ~ (dp/dH)' 1 is weakly temperature dependent. The cosine term in Eq. (0) comes from 



the invariance of the system to an overall shift in the labeling index n of the lines. It picks 
out the discrete nature of the flux lines and is crucial to the formation of a glass phase. 

Upon renormalization, one generates a term of the form V'(x, z)d z <ft, which randomly 
biases the local tilt of the flux lines. It is found that the variances of V and V are renor- 
malized in the same way, so that the inherent spatial anisotropy and frustration present 
in Eq. (p]) disappear at large length scales. It is then more convenient to work with the 
isotropic Hamiltonian, 

m= l {f (V0)2 " ^ ' W ~~ W(0(r) ' r) } ' (2) 
where r = {x, z}, /x(r) = V (r)x+V (r)z is gaussian distributed with mean zero and variance 
of the component fa 



fair)^') = a5 i:j 5 (r - r'), (3) 

and W(<j),r) oc cos[0(r) — /3(r)] is a random potential, describing the effect of a random 
phase /3(r), with 



W{(p,r)W((j)',r') = 2c/cos[0-0']5 2 (r-r'). (4) 

Denoting the bare parameters by the subscript 0, we have go ~ cr p 2 . A renormalization- 
group analysis fl4]j2lH yields the recursion relations under a change of scale by b = e l , 



dn _ 
dl ~ 

da__ 
dl ' 



0, (5) 
Ag\ (6) 



f = *9 - Cg 2 - (7) 

The coefficients A and C are cutoff dependent, however the ratio A/(kC) 2 = 2ir + 0(e) is 
universal. Eqs. (^) and (0) are valid to leading order in g, which will be sufficient provided 

T 

(8) 



2lTK 

is small. This is true even though a can flow to large values by Eq. (|6|), because the random 
potential fj, in Eq. (0) can be shifted away by the transformation 0'(r) = 0(r) — u(r), with 



kV 2 m = V /x, regardless of the magnitude of /x. The resulting potential W'(4>', r) = W((f>' + 
u,r) has the same statistics as W(<fr,r) since the latter is uncorrelated in r. Consequently, 
the flow of g cannot be affected by fi. A similar use of the statistical symmetry of 7i shows 
that the result that k is unrenormalized in Eq. (|5]) is exact |21| , |22 |. 



Clearly, e = is a special point; it defines a critical temperature T g = 4ttk, through 
Eq. (H). For T > T g where e < 0, g renormalizes to zero and at long scales the system is 
described by the gaussian part of the Hamiltonian (0) with a finite renormalized a. This 
is the flux liquid phase with the disorder causing only short wavelength modifications of 
the pure system |19| . But for T < T g where e > 0, there is a nontrivial phase controlled 
by a fixed line g*{T), with a renormalizing as in Eq. (|6|). Close to the transition, we have 
g* ~ e/C oc (Tg — T), and on scale L, a{L) « A(g*) 2 log pL. This is a vortex glass phase |23 



The existence of a perturbatively accessible fixed line allows us to study the vortex 
glass phase quantitatively. The nonrenormalization of k implies a simple form for the mean 



square thermal fluctuations of 0, i.e., ([0(r) — </>(r')] 2 ) c = T/(2ttk) log[p|r — r'|] the same 
as in the absence of randomness [J22|. The glass phase is instead distinguished by more 



strongly divergent static distortions. For example, the mean square (thermally averaged) 



displacement is (0(r) — 0(r')) m e 2 log [p\r — r'|] due to the logarithmic divergence of a [|2T] . 
However, this is not a unique feature of a glass phase, as systems with long-range correlated 
/j,'s can also give rise to anomalous mean square displacement even if g — 0, in which case 
the system is harmonic and trivial. 

We therefore consider other quantities whose behavior is unique to a glass phase. We 
first study the magnetic response of the flux array. We change the applied external field by 
an amount 5H = 5H x 5i + 5H z z, which tends to compress and/or rotate the flux array. For 
an isotropic system, the Hamiltonian becomes 

HM=H[<I>]- f h-V<}>, (9) 

where h = (5H z St + 5H x z)§ /(8tt 2 ), $o being the magnetic flux quantum. The change in 
the flux density is (d x (f>) / (2tt) , and in the "tilt angle" is (d z (p) /(27rp). The linear response 



on which we focus is Xi,j = ^(djfi)- F° r the isotropic system (H), we have Xij — X$iji an d 
the magnetic permeability is just [$q/(167t 3 )] \. Simple rescaling yields a similar result for 
the anisotropic system. 

Consider the high temperature phase where discreteness is irrelevant, i.e., g = 0. Then 
the last term in ([]) can be simply shifted away by the transformation 

0'(r) = 0(r) - h • t/k, (10) 

yielding a free energy 

F(h) = -f L i-U h- M , (11) 
and hence a response x — V K - Since the random part of F(h) is linear in h, the response 



will be sample independent as in a pure system, with (Ax) n = for n > 2 where Ax = X~X- 
This is solely a consequence of the quadratic nature of the Hamiltonian with g = 0. 

The magnetic response in the low temperature phase, where the random phase term 
in (Q) is relevant, is much more interesting. Since the transformation (pi]) does not change 



the statistics of the Hamiltonian 7i pl| , p2]| , except for generating an extra quadratic term as 



in Eq. (0), the quenched-averaged free energy Fih) is the same as for g = 0. Thus, x = V K 



independent of g. Furthermore, the average of higher order nonlinear susceptibilities (A%)' 



all vanish due to the statistical symmetry |22] . Thus average response functions are identical 
in the glass and liquid phases. This result has led some to doubt mistakely even the existence 
of the glass phase |12[ . However we will show that the glassy effects are manifested in sample- 
to-sample variations of the susceptibility and its extreme sensitivity to small perturbations. 

Let us compute the effect of the random potential W on the susceptibility variation, 
perturbatively at first. After the transformation (|10|) , the correlations between the free 
energy at two different fields hi and h 2 can be calculated to the lowest order in g. For 



AF(h) = F(h) - F(h), we have 

AF(hx)AF(h 2 ) = 2g(pL)-^ J cos [(hi - h 2 ) ■ t/k] (12) 



for a system of size L x L, with the (pL)- T ^ 2n ^ factor arising from averaging over thermal 
fluctuations. Differentiating with respect to h x and h 2 leads to nontrivial sample-to-sample 
variations of the magnetic susceptibility, with variance 

(^F=^-(p^) £ (13) 

to first order in g, with D being a sample-geometry dependent coefficient. 

For T > T g (i.e., e < 0), Eq. fll^D gives the form of approach to the asymptotic liquid 



phase where (Ax) 2 = as discussed above. For T < T g , (Ax) 2 diverges since e > 0, 
indicating the failure of the small-g expansion. Eq. ( |13|) does, however, suggest the form of 
the correct behavior: the term g(pL) e should just be replaced by the renormalized gn(L). 
Explicit computation shows that this is indeed the case. For large systems in the vortex 
glass phase, gn(L) — > g* hence we obtain a fractional variance 



^*^ = & d4) 

with D a universal geometry and boundary condition dependent coefficient, that is inde- 
pendent of nonmeasurable bare parameters such as go . For an isotropic square sample with 
periodic boundary condition, we have D rs 8ir/5. The large sample-to-sample variations 



of x indicate that the vortex glass phase is radically different from the fluid phase p4 
The size-independent variations of x i n the vortex glass phase are reminiscent of "universal 
conductance fluctuations" in disordered metals | 25[ 



Experimentally, variations of x may be obtained by measuring the magnetic response of 
one sample at different applied fields H. It will be particularly convenient to keep |H| and 



T fixed, and follow the response as the direction H is changed. The variance (Ax) 2 only 
depends on T and |H| (through k). Then as H is changed, say by rotating a sample in 
a fixed field, it effectively samples different "realizations" of the random potential, drawn 
from the same distribution since systems with different field directions H are statistically 
equivalent |]2£|| . For a system of size L x L, the free energies and hence the susceptibilities 
become uncorrelated if H is changed by an angle much greater than (pL)^ 1 as can be 
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guessed from Eq. dl2|) with k ~ §oH/p. In the glass phase, we thus expect to obtain a 
wildly varying susceptibility whose precise form is a property of the specific sample, 



but with universal statistics, in particular, (Ax) 2 /x 2 - Alternatively, one could monitor the 
magnetic noise as a function of H. This should exhibit universal variations like x(H), since 
the two are related by the fluctuation-dissipation relation. The susceptibility variations at a 
fixed T and |H| provide "magnetic finger prints" of the glassy flux phase. The reproducibility 
of the magnetic finger print for the same sample under identical conditions provides a probe 
of thermal equilibrium on long scales: Only samples small enough to equilibrate fully (see 
below) will show reproducible behavior. 

From the above discussion, it is evident that the equilibrium state of the flux array 
depends sensitively on small changes in the applied field. As argued in Refs. |L5| and 



on 



general grounds, a wide class of random systems can exhibit such sensitivity to small changes 
of a variety of parameters such as a field or temperature. Large variations resulting from 



small changes in the random potential V(r) have been studied numerically by Zhang |TJ 
for a single flux line. In the remainder of this paper, we analyze explicitly the effect of such 
a small change in the random potential for the 1+1 dimensional flux array. Sensitivity of 
the array to small temperature changes can be analyzed similarly. We merely quote the 
analogous result for this somewhat more complicated case. 

We consider two noninteracting flux arrays, 0(r) and 0(r), in two different realizations 
of the random potential, {/j,(r),W(<p,r)} and {jS(r), W(<j>, r)} respectively. We take the 
random potentials to be statistically equivalent but slightly different from each other, so 



that /x(r)/x(r') is given by Eq. (|3|) and W((f>,r)W((p' ,r') given by Eq. (Q). However the 
cross-correlators are 



/iitr)/!^!-') = a8ij8 2 (r - r' ), (15) 
W{4),r)W{ft,Y') = 2£cos[0-0 , ]5 2 (r-r'), (16) 

with the bare values So < &o and go < go- The renormalization group recursion relations 
Eqs. (|) - (0) must be unchanged as the systems are uncoupled. However, the cross corre- 



lations renormalize as 

§= ^\ (17) 
§=eg-Cgg, (18) 

with e = e+(a-a) /(2ttk 2 ) . 

To investigate the effect of weak de-correlation of the random potentials in the glass 
phase, we linearize the recursion relations around the vortex glass fixed point g*. For small 
de-correlation S a = a — a <C 1 and S g = g — g <C 1, we have 

<f= 2Ag'S g , (19) 

*= SS?^ (20) 

This flow has one positive eigenvalue 

X 5 « g*yjA/7TK 2 = V2e. (21) 

Therefore, infinitesimally small de-correlations in the bare random potential grow under 
renormalization. On long scales L ^> L5 ~ 5 _1 / A *, with 5 a linear combination of 5 a and 8 g , 
g(L) vanishes and d{L) saturates. The two systems then appear substantially different and 
will have essentially independent susceptibilities, with Ax{L)Ax{L) — > for large L. There 
will, however, be residual cross correlations associated with the finite renormalization of a. 

These effects can best be probed by changing the temperature of one sample slightly 
by ST. The same exponent A5 in Eq. (|21|) controls the crossover, and for system sizes 
L > L 5 ~ (5T)- l ^ s , x(T) and x{T + ST) will be roughly independent. If ST < T 9 - T, 
the temperature dependence of x w ih probe statistically similar variations of x as did the 
field direction dependence of x(H). 

Physically the source of the sensitivity to H and T changes are quite different. The 
former is due to the changes in mean position Of the lines while the latter is more subtle: 
It is caused by the entropic contributions to the free energy, which drastically changes the 
effective random potential on long scales. Although this has been predicted for a variety of 

8 



random systems [ I5|,18| and supported by numerical and approximate renormalization group 
calculations, this to our knowledge, is the first time an analytic calculation has yielded the 
hypersensitivity to temperature changes. 

We close with a comment on dynamics: Recently, a number of authors have claimed 
that free energy barriers in this system grow as various powers of logL |7|||,[nj. An explicit 
dynamic renormalizat ion-group calculation fli~l~l , |27| found that the dynamic exponent z ~ 2 + 
1.8e for T < T g , yielding a nonlinear resistivity, dS/dJ ~ J a9e , where £ is the EMF generated 
by a uniform current J applied perpendicular to the (x, ,z)-plane. However, because the 1+1 
dimensional vortex glass phase is controlled by a finite temperature fixed line rather than a 
zero temperature fixed point, the barriers are not well defined by the form S(J) found. In 



the limit T — > 0, however, one finds |28] z ~ which can then be correctly interpreted 
as barriers growing as logL. 

The dynamics can also be used to probe the length dependence of x- At finite frequency, 
to, scales of size l u ~ oj~ 1 I z are probed. Since the susceptibility x( u ) f° r eacn correlation 
volume 1^ will be essentially independent, the variations in x(u;) for a sample of size L x L 
will be Ax(cj) ~ l u /L ~ cu^^/L, crossing over to the static result only when ~ L. 

In this paper, we have analyzed some of the glassy properties of randomly pinned flux 
arrays confined to a plane. In the vortex glass phase, the magnetic susceptibility is found 
to be strongly dependent on the external field, temperature and specific sample, exhibiting 
variations with universal statistics. 
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